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The recent excess in the CMS measurements of eejj and evjj channels and the emergence of 
PeV cosmic neutrino events at the IceCube experiment share an intriguing implication for a lepto¬ 
quark with a 600-650 GeV mass. We investigate the CMS constraints on the flavor structure of a 
scenario with the minimal leptoquark Yukawa couplings and correlate such a scenario to the reso¬ 
nant enhancement in the very high energy shower event rates at the IceCube. We find for a single 
leptoquark, the CMS signals require large couplings to the third generation leptons. This leads to 
an enhancement in the ^v-nucleon scattering cross-section and subsequently more v T events at PeV 
energies. However, a visible enhancement above the Standard Model scattering would require a 
leptoquark Yukawa coupling larger than one that can be easily tested at the upcoming LHC runs. 


I. INTRODUCTION 

CMS has measured p] a 2.4tr excess in eejj and evjj 
channels, compared to the earlier null results from the 
CMS !2] and ATLAS [3] Collaborations. The CMS ex¬ 
cess has indicated a best fit for a leptoquark with 650 
GeV mass pQ. This interesting possibility has been in¬ 
vestigated in recent works [3j. Through its recent ob¬ 
servation of high energy neutrinos 0 0, the IceCube 
has potential to test leptoquark models. Interestingly, 
at the PeV energy range, the neutrino flux (up to un¬ 
certainty due to the limited statistics) may be fit within 
the picture of a leptoquark resonance [7] . Additional ex¬ 
planations to the IceCube neutrino flux include a new or 
modified cosmic source flux [3J, new neutrino flavor(s) 0, 
colored-neutrino resonances [ini, a Glashow resonance at 
the next-to-leading order Standard Model explana¬ 
tions [12], as well as other possibilities [13] . 

In the Grand Unified Theories (GUTs) [H| such as 
SU( 5), S'O(IO), and E e , etc, the lepton (L) and baryon 
( B ) number symmetries are violated explicitly, although 
the B — L symmetry can still be preserved. Thus, lep- 
toquarks generically exist which mediate transitions be¬ 
tween quarks and leptons. In addition, these leptoquarks 
will not contribute to gauge anomalies since they are 
scalars. So they provide simple extensions to the Stan¬ 
dard Model (SM). 

This paper assumes a single scalar leptoquark that 
explains both signals from the IceCube and the CMS, 
and shows that these current experimental hints demand 
a specific flavor structure in the leptoquark hypothesis, 
which can make predictions for the future PeV events at 
the IceCube and will be probed efficiently by the upcom¬ 
ing LHC data. 

Section [TT] gives the minimal flavor assignment for 
the leptoquark couplings and possible constraints. Sec¬ 
tion |III| and |IV| present the leptoquark’s signal rates at 


the CMS and the IceCube. We then conclude in Sec¬ 
tion El 

II. A MINIMAL LEPTOQUARK SCENARIO 

We follow the notations of Ref. m and focus on a 
single scalar ‘SI’ type leptoquark S whose SU(3)c x 
SU(2 )l x U(1)y quantum number is (3,1,1/3). Its 
Yukawa couplings to the SM quarks and leptons are 

c => yijQi,L iT 2ij,LS + y*ul R ij,R.S + 

+y'ij DF 0 C i,R d j,RS^ + c.c. , (1) 

where the subscript i,j denote the flavor generations of 
quarks and leptons. The matrix 1 T 2 swaps the two isospin 
elements in Ip. To solve the proton decay problem, con¬ 
sider the Z 3 symmetry as follows 

Qi,l —■► u 2 qi , Ui }R —► uj 2 Ui tR , d itR —► u) 2 di, R 
U. -A u)li , ei :R -A we^ij , (2) 

where to 3 = 1. Thus, we forbid the diquark couplings 
yf® and y R j D , and then escape the constraints from the 
proton decay experiments. 

While a leptoquark can generally couple to all three 
generations of quark and lepton flavors, a minimal cou¬ 
pling scenario can be favored: 

(i) The IceCube neutrino events occurs via neutrino 
scattering off the nucleons in ice, which demands left- 
handed coupling to the valence (first generation) quarks. 

(ii) Due to the high energy nature of incoming cosmic 
ray neutrinos, an s-channel leptoquark resonance cross- 
section is suppressed due to the contribution to the lep¬ 
toquark decay width from any couplings other than that 
in (i), i.e., all right-handed couplings y R and any left- 
handed y L to the second or third generation quarks. 


2 


(iii) The lack of muon tracks among the highest energy 
events disfavors its couplings to y and v^. 

(iv) The leptoquark Yukawa couplings to heavy quarks 
introduce constraints arising from l —> I'j decays. 

Keeping the conditions (i)- (iv) in mind, we require 
our minimal flavor structure to consist of nonzero y R 3 
and y R i- The leptoquark interaction Lagrangian is 

£lq = 2/i3<2i, L^h, lS + y R ^ R e R S 

+ U^rTrS + 0(e) + c.c. , (3) 

where <71 = (u,d),l 1 = {v e ,e) and l 3 = (v Ti t). y 33 is re¬ 
sponsible for the ^ T -nucleon collision while 2 / 21 , 2/23 m ust 
be included to balance the CMS signal rates in the eejj 
and evjj channels. y R does not mediate neutrino-nucleon 
scattering. We choose the right-handed charm quark in 
stead of up quark so that the right-handed couplings do 
not cause a large enhancement to the rare 7r + — > e + v de¬ 
cay m and a large radiative correction to the up quark 
mass, when the left-handed coupling to the 1st gener¬ 
ation quark doublet is also present. The couplings in 
Eq. [3] can lead to D + —>■ e + v decay, yet given the uncer- 
tianty in the D + meson form factor QjS], it can tolerate 
relatively large couplings y ~ 1. O(e) denotes for other 
terms with couplings that are negligible in the quadrature 
sum J2ij(2\y L \ij + \y R \ij)- To generate an eejj signal, in 
principle there can also be y[\ . However y 31 is also not 
as effective as y R \ to explain the CMS signal ratio which 
will be shown shortly. 


III. THE COLLIDER SIGNALS 


The CMS experiment has observed [Tj 36 eejj and 18 
evjj events at 19.6/6” 1 whereas 20.49 ± 2.14(sfaf) ± 
2A5(sys) and 7.54 ± 1.20 (stat) ± 1.07(s2/s) events are 
expected from the standard model (SM) background for 
these final states respectively. A best-fit with 650 GeV 
leptoquark would require about 20% decay branching ra¬ 
tio (BR) in the S —> ej channel, subject to possible vari¬ 
ation in the leptoquark pair production cross-section at 
the next-to-leading order [Tj 1 (NLO) and the experimen¬ 
tal event selection cut efficiencies. For slightly lighter 
leptoquarks, the BR ej - can be lower due to a greater pro¬ 
duction rate. 

Converting the CMS event rates into signal cross- 
sections, we get 


LQ 

eejj 

LQ _ 


0.79 ± 0.33 (stat) ± 0A3(sys) fb , 
0.53 ± 0.22(sfat) ± 0.055(s2/s) fb , 


(4) 


1 BR; 


■ej 


*eejj 


'eejj _ _ 

^ei/jj 2 BR^BR e j .'l e.wj j 


= 1.5 ±0.9 , 


1 As a note, in case of a very large leptoquark couplings to the 
1st generation quarks, t-channel lepton diagrams should also be 
included beside QCD processes. 


where we have added the subdominant systematics in 
quadrature for r ev . A denotes the cut efficiency in 
each channel. We find both A eej j and A ev jj ~ 60% 
using MadGraph v5 and electron/jet reconstruction by 
Pythia+PGSj for the selection cuts listed in Ref. [T]. r 
leptons only have a few percent chance to fake an elec¬ 
tron event and can be neglected. Ignoring the quark and 
lepton masses, the leptoquark decay width is 

r = ^yE 2 %) 2 + (»S) 2 . W 

i,3 

where the factor of 2 arises from the left-handed coupling 
giving rise to both S -A Iq and S —> viq’ decays. With a 
minimal set of { 2 / 13 ,2/iii 2 / 21 } couplings, we have 

_ 1 „ (2/i L i) 2 ±(2/H) 2 
el/ 2 (2/u) 2 + (2/13 ) 21 

RR = (2/n) 2 + ( 2 / 21 ) 2 

2(2/^ + 2{y R 3 Y + (y*y + {y*Y 

Both 2 /n and 2/21 terms contributes to BR ej . Since 
the CMS has observed less signal evjj than evjj , from 
Eq. ([t]) it is clear that the 2/21 term is more favored com¬ 
pared to the 2 /n term which makes it difficult to keep 
BR vj down. We, therefore, drop y 31 from now on. Al¬ 
though the central value prefers a larger y R , due to the 
still significant statistical uncertainty, 2/21 ~ 2/13 is still 
allowed, and two benchmark points are given in Table |T] 


AIlq 

I 2/2l/2/l3 I"" 2/23/2/13 

BRej 

BR„j 

BR rj 

600 

1 

2 

20% 

20% 

60% 

650 

0.57 

0.29 

20% 

35% 

45% 


( 6 ) 

(7) 


TABLE I: The benchmark points for the leptoquark Yukawa 
coupling assignment schemes that are consistent with the 
CMS measurements. 

The scalar leptoquark pair production is QCD domi¬ 
nant, and its next to leading order (NLO) cross-section is 
given in Ref. E 3 ■ At the benchmark points, a 600(650) 
GeV leptoquark has a cross-section of <r 0 = 20(11) fb. 
With a 60% selection cut efficiency, it leads to the after¬ 
cut signal (Teejj = 0.5(0.3) fb and (J e vjj = 0.96(0.9) fb, 
within 2a uncertainty of the CMS observed event rates. 

It is worthwhile to notice that a sizable 2/23 is needed 
to prevent y 33 from dominating all of the non-ej decay 
branching ratio, in which case a large signal a ev jj may 
easily be ruled out by the current CMS measurement. 


IV. ICECUBE NEUTRINO EVENTS 

At the IceCube detector, the incoming neutrinos with 
PeV energy would have enough center-of-mass (COM) 
energy ^ /2m p E v > Mlq to trigger a resonant s-channel 
leptoquark exchange. Due to the threshold neutrino en¬ 
ergy for the leptoquark resonance, the neutrino-nucleon 
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scattering cross-section deviates from that in the SM (as 
in Fig. [Tj) , and can generally explain the emergence of 
the few PeV events with a coupling y L ~ 1 |7j. In com¬ 
parison, while an explanation using the SM only can also 
account for the observed number of PeV events due to 
its power-law spectrum which follows the shape of the 
incoming cosmic neutrinos, it would be difficult to de¬ 
velop features at a particular energy scale in the event 
energy spectrum. Admittedly the current data is statis¬ 
tics limited and the ‘discontinuous’ spectral feature needs 
future measurements, the leptoquark resonance is a very 
interesting interpretation, and can be predictive when it 
is used to explain the CMS excess as well. 

If the leptoquark interactions only involve light quarks, 
the differential resonance cross-sections are simplified 
compared to that in Ref. [TSj, 


dcr NC / cc 

dy 

n cc 

n vc 


* D{Mj Q /s) 
2 K s 

ylivl + VR ) 

2 y | + v 2 r ’ 
Vl 

2 Vl + Vr ’ 


( 8 ) 

( 9 ) 

( 10 ) 


where y = Ei/E v is the daughter lepton energy frac¬ 
tion. D is the averaged 2 nucleon’s parton probabil¬ 
ity distribution at the resonance COM energy fraction 
s/s = M\q/s. Depending on the leptoquark’s decay 
channel, the analogue to the SM’s charged-current (CC) 
final states S —>• Iq deposit all the incoming energy into 
visible showers; in a neutral current (NC) S —> uiq' event 
the invisible neutrino hides away a fraction of the total 
energy. Therefore, the CC events are better candidates 
for the several highest-energy IceCube shower events, 
while the lower, sub-PeV ones could be a continuum due 
to NC events. 

In terms of the couplings shown in Table [TJ we have 

Vl = (2/13 ) 2 and Vr = (vi i ) 2 + (z/H) 2 - :t is clear that the 
CC process is less suppressed from sizable right-handed 
couplings compared to the NC process, as S —>■ viq' 
is strictly left-handed. Fig. [l] shows the SM and the 
leptoquark-enhanced neutrino cross-sections. Due to the 
appearance of the leptoquark resonance, the total cross- 
section exhibits a rise near the threshold neutrino energy 
E th = Ml Q /(2M p ).. 

The cosmic neutrino flux is assumed to be isotropic, 


dcj) 

dE v 


= 3 </>o/i 


E v \~ 7 
100 TeV ) 


( 11 ) 


where /,; is the fraction of neutrinos of the ith flavor. The 
observed spectral parameters are given as 7 = 2.6 ± 0.15 
and <j >0 = (2.3 ± 0.4) x 10 -18 GeV _ 1 s _ 1 cm _ 2 sr _1 in 


2 With a 5:4 proton to neutron ratio in ice. The d, quark’s proton 
distribution is taken to be same as the u quark inside a neutron. 



FIG. 1: (Left) v T neutrino-nucleon scattering cross-section 
from the SM and a 600 GeV leptoquark resonance with y± 3 = 
1. Both NC and CC interactions are included. v £ and v p 
do not receive enhancement under our flavor assignment, and 
the leptoquark contribution to v T is subdominant. 


Ref. IS], which also reported good consistency with the 
conventional neutrino flavor ratio (| : | : |), and a best, 
however unphysical, fit at (0 : 0.2 : 0 . 8 ) with the cur¬ 
rent data. While muon track events only originate from 
ry,, the shower-events include contributions from the NC 
scattering from all favors as well as the CC scattering 
from both v F and v T . For shower events, we only con¬ 
sider the ‘contained’-type of events in which the shower 
develops inside the Icecube’s detector array volume. The 
contained shower event spectrum is 


dN_ 

~dE 


= Nt 


targets 


poo 


da NC (E v ) difiv 


+ 


I £ dE dE u 

da cc (E„ r ) dcp Vr 


dE„ 


IE 

cc. 


+ a cu (E) 


dE 

~dE 


dE , 


dE„ T 


V T 


( 12 ) 


The leptoquark component in the NC scattering cross- 
section can be easily derived from Eq. (10) as dcr/dE = 
1 /E v ■ da(E u )/dy, where E is the visible energy deposit. 
The v e CC scattering deposits 100% of the neutrino en¬ 
ergy into showers, while for v T CC scattering, a loss in 
the energy deposition occurs due to the secondary neutri¬ 
nos that emerge from r decays. The CC showers from v e 
and v T have unequal hadronic contents and consequently 
generate different amount of light yields. A full shower 
acceptance simulation is beyond the scope of this paper, 
and we adopt the approximation given in Ref. m that 
100% of v e CC shower is measured, while for v T it is ap¬ 
proximated as E ss 0.5E„ t , or da/dE = a cc (E V t )- 5(E— 
E Vt / 2). Similarly, for the mainly hadronic NC showers, 
we use da/dE = a NC (E„) ■ 8{E — 0.2 E„) for all flavors. 
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Just like the SM, both CC and NC processes con¬ 
tribute. However, unlike the SM, the leptoquark res¬ 
onance only occurs involving a neutrino and a quark, 
while the CP-conjugate process, the anti-neutrino res¬ 
onance scattering is highly suppressed due to the anti¬ 
quarks being only sea quarks inside the nucleons. The 
parametrizations given in Table [I] only needs to include 
the leptoquark enhancement in the v T scattering. 



E [PeV] 



E [PeV] 


FIG. 2: The shower event rates in the southern hemisphere. 
The leptoquark contribution assumes Mlq = 600 GeV and 
2 /f 3 = 1.4. The first panel assumes a (| : | : |) flavor ratio 
while the second panel assumes the best-fit (0 : 0.2 : 0.8) 
ratio, which makes the LQ signal more pronounced with less 
SM CC showers. 


Fig- i shows the leptoquark’s enhancement to PeV 
shower events in the southern hemisphere. We show both 
the conventional ( 5 : 5 : 5 ) and the best-fit ( 0 : 0 . 2 : 0 . 8 ) 
flavor ratio. In the latter case the leptoquark signal be¬ 
comes enhanced due to its mostly third flavor coupling. 
Since the leptoquark resonance only occurs in the PeV 
events, the energy dependence in the effective volume 
is dropped, which would be more significant around or 
below tens of TeV. A constant (saturated) detector ac¬ 
ceptance is assumed and the SM event rate is calibrated 
to the sub-PeV data. The red shaded region assumes 
a coupling size ?/f 3 = 1.4 for the first benchmark sce¬ 
nario in Table [TJ The second benchmark results in a very 
similar signal and is not shown. The majority of the lep¬ 
toquark signals occur via u T CC scattering, which would 
indicate itself through a rise in the fraction of rs in the 
PeV events. 


Interestingly, the significance of the leptoquark contri¬ 
bution can be sensitive to the neutrino flavor composition 


under the fractional energy deposition scheme in Eq. 12 


where v e CC presumably deposits 100% neutrino energy 
and becomes the most efficient contributor to the PeV 
event rate because it picks up cosmic flux at a lower en¬ 
ergy, which falls steeply as an E~ 26 power-law. As a 
result, the SM shower rate becomes dependent on the 
v e composition. As shown in the lower panel in Fig. [2} 
the (0 : 0.2 : 0 . 8 ) flavor ratio, which would require be¬ 
yond SM sources, leads to a smaller SM rate relative 
to an enhanced (by a higher v T composition) leptoquark 
rate, making the leptoquark more visible. The lower v e 
CC contribution also causes a difference in the effective 
exposure when calibrating to the sub-PeV SM events, 
which is 650 Megaton of water-equivalent times year for 
the ( 5 : 5 : 5 ) flavor ratio and 1.3xl0 3 Mton-year for 
(0 : 0.2 : 0 . 8 ), at the central values of parameters used in 

Eq-EU 

As a caveat, this effect from neutrino composition 
might be over-estimated because v T CC and the NC 
events in principle can deposit more than their average 
fractions of energy in dcr(E v )/dE. A proper detector 
simulation would then be required, especially in the case 
of I'e-depleted sources, to make a more precise predic¬ 
tion. Also, while a leptoquark enhancement to v T scat¬ 
tering may make the source appear to have more v T in 
a SM-only calculation, the relatively small leptoquark 
event rate (as seen in the upper panel of Fig. [2]) with a 
(| : | : |) flavor is insufficient to boost the cosmic source 
into a (0 : 0.2 : 0.8) appearance. So we utilize the latter 
only to demonstrate the leptoquark’s sensitivity on the 
flavor composition of the neutrino source flux. 

In Fig. [$J we illustrated with |yf 3 | 2 = 2 . The lepto¬ 
quark event rate at other coupling values can be simply 
scaled by \y\ 2 . The enhancement raises the event rate 
over a wide range of energy. More events are expected 
just above the threshold, as can be seen in Fig. [T] In 
the high energy limit, the leptoquark contribution ap¬ 
proaches to a constant enhancement and the total event 
spectrum would return to a power law which follows the 
shape of the incoming cosmic flux. 


V. DISCUSSIONS 

It is clear that the presence of a leptoquark helps to fit 
the IceCube data even though the inclusion of the lepto¬ 
quark alone may not fully explain the shape of the PeV 
shower spectrum with only a few events at present, where 
the SM still plays a major role. If excessive PeV shower 
events keep occuring in the future data, a change in the 
power-index of the event energy spectrum at the PeV 
energy scale can be interpreted as a leptoquark signal. 
IceCube’s discovery analysis uses a very conservative se¬ 
lection criteria to extract a very low backgrounds sample, 
given that the flux was unknown. With a better under¬ 
standing of the signal, re-optimized selection criteria can 
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improve event statistics significantly, for example, by re¬ 
laxing the containment criteria. Further, the interest in 
gigaton size detectors has substantially increased, with 
several projects targeting the high energy region, e.g., 
IceCube-Gen2 [20], KM3Net in the Mediterranean 2T| . 
and Gigaton Volume Detector (GVD) [22] in lake Baikal. 
By the end of this decade, IceCube alone will triple the 
sample size (by the most conservative estimates) in the 
PeV region due to the increased exposure. In the early 
next decade, using upgrades like Gen2 alone or in combi¬ 
nation with GVD and KM3Net, one could easily obtain 
a ten fold increase in the annual event rate and a pre¬ 
cise measurement of the PeV spectrum would then allow 
us to distinguish between conventional and leptoquark 
scenarios. 

If the PeV events were enhanced by a leptoquark 
around 600-650 GeV, the CMS data would require this 
leptoquark left-handedly couple to the 3rd generation 
leptons and the first generation quarks, also with sizeable 
right-handed couplings, as represented by the benchmark 
points. This would predict more v T events at PeV. How¬ 
ever, a fairly large coupling yi, > 1 is required for the 
leptoquark contribution to vN scattering to be compa¬ 
rable to that from the SM. The upcoming LHC data in 


the eejj and evjj channels can easily verify or exclude 
leptoquark couplings of such a magnitude, as well as the 
flavor structures shown in Table |H 
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